Model for Density Waves in Gravity-Driven Granular Flow in Narrow Pipes 
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A gravity-driven flow of grains through a narrow pipe in vacuum is studied by means of a one- 
dimensional model with two coefficients of restitution. Numerical simulations show clearly how 
density waves form when a strikingly simple criterion is fulfilled: that dissipation due to collisions 
between the grains and the walls of the pipe is greater per collision than that which stems from 
collisions between particles. Counterintuitively, the highest flow rate is observed when the number 
of grains per density wave grows large. We find strong indication that the number of grains per 
density wave always approaches a constant as the particle number tends to infinity, and that collapse 
to a single wave, which was often observed also in previous simulations, occurs because the number 
of grains is insufficient for multiple wave formation. 

PACS numbers: 45.70.Mg, 83.10.Rs, 47.11. Mn, 47.57.Gc 
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I. INTRODUCTION 

Transport of dry granular media through pipes and 
channels is a problem of fundamental interest which is of 
considerable industrial importance pLj], and it has been 
studied intensely both experimentally and theoretically 
in recent decades . Granular media behave radically 
different from both liquids and solids Q and in granular 
pipe flow driven either by gravity or pressurized gas, non- 
linear dynamical phenomena such as clogging (density 
waves) are observed, but not yet well understood. 

We consider dry grains falling inside a pipe. In this 
system, encountered industrially, e.g., in emptying of si- 
los and transportation of sand or powder, there are three 
main mechanisms of interaction: (a) collisions between 
grains, (b) collisions between grains and walls, and (c) in- 
teraction of the grains with air in the system. Assuming 
the pipe is narrow, we approach this problem by means 
of a simple one-dimensional model with periodic bound- 
ary conditions in which collisions are modelled by means 
of two coefficients of restitution fi and v, corresponding 
to the collisions of mechanisms (a) and (b), respectively. 
We demonstrate that this is sufficient to observe the for- 
mation of density waves. 

Granular media has frequently been investigated by 
means of numerical simulations. Features such as clus- 
tering through dissipative collisions [1, 0] and inelastic 
collapse d, 1^ are among those reported. Various driv- 
ing mechanisms have been considered, which can roughly 
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be divided into two categories: vibrating walls 
and interior heating |13l - [l8| . The case of an emptying 
hopper has been studied, see e.g. 0, [13], as has the de- 
compaction transient associated with the release of an 
initially dense packing [2l| . 

The particular case of gravity-driven granular flow 
through channels and pipes has been studied for some 
time, both experimentally and theoretically. In previ- 
ous experimental investigations [22| - [28j the importance 



of the presence of air in the flow has been emphasized, 
and the formation of density waves in falling sand has 
been explained with primary reference to the air-grain 
interaction. Although clearly of importance in real sys- 
tems, we demonstrate that the introduction of air in the 
model is not necessary to observe density waves. Our 
simulation suggests that such flocculent behaviour, while 
certainly influenced by the presence of air, can occur also 
in evacuated pipes, provided dissipation from collisions 
between grains and walls is faster than dissipation from 
grain-grain collisions. 

Previous simulations of granular pipe flow have also re- 
ported density waves [3, [29l - [33 | . Several mechanisms for 
velocity dissipation through inelastic collisions, damping 
and static friction have been employed, see, e.g., |3d|. 
Early papers [l,|2l| were unable to study sufficiently large 
numbers of grains and collsions to reach a state indepen- 
dent of initial conditions. The model used by Lee [s^] 
was a 2-|-l-dimensional time-driven molecular dynamics 
(MD) simulation, with at least eight significant param- 
eters. In the model of Peng and Herrmann [31|, |3^, a 
2+1 dimensional lattice gas automaton is used in which 
various events are assigned probalistic collision rules. An- 
other model was proposed by Liss, Conway and Glasser 
[ssj . Again, this is a 2+1 dimensional model with more 
complicated collision rules than ours. In their simulation 
the grain-wall coefficient of restitution was found to be 
of little importance. In all these models stable density 
waves were seen. 

Our model is one-dimensional and much simpler, yet 
reproduces similar qualitative phenomena, dictated es- 
sentially by two coefficients of restitution whose inter- 
pretation is physically transparent. 

We present our model and the results of our simula- 
tions and demonstrate that the model has two clearly 
distinct regimes in the parameter plane, one in which 
density waves form (flocculent regime) and one where no 
such are present (gaseous regime). Average number of 




FIG. 1: Spatiotcmporal plots of grain trajectories, (a) Gaseous regime /i — 0.8, = 0.3. (b) Split velocity transition point 
fi — = 0.55. (c) and (d): Different multi-wave patterns, (/i, i^) = (0.55,0.65) and (0.30,0.85), respectively, (e) and (f): 
Collapse to a single wave; = (0.70,0.95) and (0.15,0.85), respectively. Further details in text. 



grains per density wave and flow rate as functions of tlie 
two coefficients of restitution are studied. 



II. THE MODEL 

The pipe of our model is sufficiently narrow to pre- 
vent grains from passing each other, and sufficiently wide 
to allow them to fall approximately freely between col- 
lisions. As a consequence, grain-wall collisions are trig- 
gered by, and occur immediately subsequent to, grain- 
grain collisions. The pipe is therefore assumed to be too 
narrow for arching effects to play a role. We thus propose 
that the system is essentially one-dimensional in behavior 
and may be captured qualitatively by a one-dimensional 
model where collisions with walls are incorporated into 
the grain-grain collision rules themselves. 

The set-up is the following. Let N grains move in one 
dimension so that grain i has position Xi and velocity Vi . 
We use periodic boundary conditions, with grains falling 
beneath x = reinserted a,t x = L, where L is the length 
of the pipe. The grains are accelerated by a constant 
gravity g towards x — 0. In the dilute limit dN <C 
where d is the diameter of a single grain, we find that 
the behaviour is independent of d, and we set d = for 
simplicity. 

As mentioned, the collision rules employed are based 
on the idea that the falling grains lose energy due to two 

different types of collisions, with the walls and with each 
other. Grain-grain collisions are modelled by reducing 



the relative velocity of the colliding grains through a co- 
efficient of restitution fj, e [0,1]. Thus, v^^^^ = — M^Rab 
where VB.,a.b = Vn — Vh, grain 'a' is directly above grain 
'b', and we use superscripts ' — ' and to denote times 
just before and after a c;ollision, respectively. 

Grain-wall collisions tend to decrease the average ve- 
locity of the system, = -^"^ X^ili ^'i; which would 
otherwise diverge in time. Instead, Vav eventually reaches 
and fluctuates around a constant value at which dissi- 
pated energy is in equilibrium with energy gained from 
the gravitational field. We model the interaction with 
the walls at each collision by a coefficient of restitution 
v e [0, 1] which dissipates grain energy by reducing the 
center of mass (CM) velocity of the two colliding grains 
according to v^y^ — i^f^M ^b' where the CM velocity 
of the pair is I'cM.ab = ^(^'a + ^^b)- These rules combine 
to 

=\{v - /i)ua + \{v + M)^'b (la) 

■"t =\{^ + lA^Z + \iy - ^J)^h ■ (lb) 

Note that introducing a restitution of vcm, momentum 
is transfered for each collision from the grains to the (in- 
finitely massive) pipe, and the momentum of the grains 
alone is not conserved. This causes the flow to approach 
a constant flow rate as observed in experiments. 

While a collision changes the velocities of just two 
grains, the relative and CM velocities of three pairs of 
grains are affected: the colliding pair plus the neighbor- 
ing pairs above and below. In particular, a collision in- 
creases the relative velocity of the pair above the collision 
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in proportion to the lost CM velocity of the colliding pair. 
Thus, relative velocities are constantly regenerated, and 
the phenomenon of inelastic collapse [3, Q is avoided for 
/i > and V < 1. 



III. SIMULATION AND RESULTS. 

We use event-driven molecular dynamics for our sim- 
ulations. Each grain is given an initial velocity which 
is assigned randomly according to various initialization 
schemes. The grains are also given equidistant starting 
positions between and L. 

Simulations are run for different grain numbers N at 
constant grain density p — N/L. After a 'thermaliza- 
tion' time the system reaches a steady state where the 
flow rate of the whole system fluctuates about a mean 
value. Initializing the simulations with different initial 
conditions yields the same dynamic steady state, charac- 
terized by average flow rate, collision frequency and the 
average number of density waves. 

In Fig. [T] we plot grain trajectories in space and time. 
In these plots we have used, in arbitrary units, L = l,g = 
0.01 and N = 100. In these simulations initial velocities 
between and 1 were drawn randomly from a uniform 
distribution. 

We observe two distinct regimes in the fj,, v plane. 
When V < p, dense regions tend to spread out. In this 
'gaseous' regime, shown in Fig. [1^, no steady density 
waves are observed. Whenever v > fx, density waves are 
observed, and the transition from gaseous to flocculent 
behavior a.t fj, = v is sudden. 

For the special case 11 = v seen in Fig. [T}), the eqs. JT]) 
simplify to t;+ — jiv^ and — fiv~ . Now the veloc- 
ity of grain 'a' after a collision is only a function of the 
velocity of grain 'b' before the collision, and vice versa. 
After 'thermalization' we observe that the grain veloci- 
ties are organized about two values: a lower velocity for 
the grains whose last collision was with the grain below 
it and a higher velocity for those which last collided with 
the grain above it. 

While the qualitative behaviour in the gaseous regime, 
p > v, is insensitive to the values of /j, and i/, the situa- 
tion in the flocculent regime is quite different. As exem- 
plified for the parameter pairs in Fig. [1] the wave pat- 
terns vary strongly in this half of the parameter plane. 
Both the magnitude and velocity of the density waves 
depend sensitively on the values of /i and v. Moreover, 
the qualitative picture varies from a few large and sta- 
ble waves with approximately constant velocity to many 
small and volatile waves which emerge, merge and dis- 
solve and vary greatly in velocity even for a single set of 
parameters. Two examples are shown in Fig. [TJ; and d, 
but a multitude of different multiwave patterns may be 
observed. 

Beyond rescaling the time axis, the values of L and g do 
not affect the wave patterns, so long as the grain number 
N is large compared to the average number of grains per 
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FIG. 2: (Color online) (a) Number of grains per waves and (b) 
average grain velocity as functions of particle number for 
three sets of parameters (same in both panels). The dashed 
line in the above panel is A'^ = N^, i.e., a single density wave 
is present. 



density wave. For TV below some /z,z^-dependent thresh- 
old, the grains will gather in one or a few stable or 
metastable waves, as exemplified in panels (e) and (f) 
of Fig. n 

When N is increased beyond this threshold, however, 
we find that the number of density waves increases lin- 
early, while the number of grains per density wave remain 
constant. Average flow rate and the average velocity of 
density waves also remain independent of TV in the limit 
of large N. This represents the asymptotic limit of our 
system, where there is no dependence on density p, which 
we illustrate for three pairs of coefficients in Fig. [2j Due 
to limitations of computer time we have not investigated 
this limit when iV^, grows beyond ~ 400, which happens 
in the limit of either small /x, or w 1, but we conjec- 
ture that such an asymptotic regime always exists when 
p and v are both on the open interval (0, 1). Single- wave 
pictures such as panels (e) and (f) of Fig. [TJ seen also 
in various other simulations |30l433l |. thus appear in our 
model simply because there are not enough particles in 
the system to form more than one wave [36[. 

In Fig. [3^ we estimate the number of grains per wave 
in the flocculent regime, N^, after 'thermalization' (cal- 
culated as N divided by the number of density waves) in 
the asymptotic (large N) regime. In practice we choose 
N large enough that the observed value of iV^, is constant 
with increasing N. We have used N ranging from 500 to 
4000 in different areas of the plotted region. TV^, varies 
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FIG. 3: (Color online) (a): Plot of log A'^j^, where A'^u; is number 
of grains per density wave as a function of /i and u in the 
asymptotic limit, (b): Plot of average flow rate as a function 
of fi and v in the asymptotic limit. In both panels, encircled 
letters refer to parameter pairs in Fig. [1] 



nonmonotonously and spans several orders of magnitude. 
Indeed, when fi decreases below 0.2 or v approaches 1, 
both iV^ and 'thermalization' time diverge rapidly, mak- 
ing calculation in this region expensive. 

In Fig. [3)d the average flow rate is given for all values 
of fi and v. In the gaseous regime we see a monotonous 
variation supporting our observation that there is no sig- 
nificant parameter dependency on the structures formed 
in this region. Within the flocculent regime the pattern 
is similar to that observed for the number of grains per 
density wave. 

Counterintuitively, flow rate is highest in areas of the 
/i, v plane in which the density waves are large, and con- 
versely, smaller waves correspond to a low flow rate. 
Since wave velocity is always lower than the flow rate, 
each grain will approach a density wave from above, col- 
lide its way through it, and fall into a low-density area 
beneath it once more before meeting another wave. Al- 
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FIG. 4: Frequency spectrum of the wave pattern of Fig. [TJl. 
The graph is the average over the power spectra of 200 sub- 
samples of a long time series of « 10^ data points. 



though large waves move more slowly than smaller waves, 
we find that this is more than compensated by the pres- 
ence of large low density areas between waves in which 
the grains can fall freely. Hence the total fiow rate is 
governed by the length of these acceleration stretches 
whereas we find no clear connection between flow rate 
and wave velocity. 

We have analyzed the power spectrum of a typical wave 
pattern (that of Fig. [TJi) shown in Fig. HI The analysis 
was performed by recording density in a region of length 
L/10 as a function of time. With sampling frequency 
10 per time unit, about 10^ data points were recorded. 
The series was then divided into subsamples which were 
Fourier transformed individually, and the average of the 
power for all sub-samples were then taken (different num- 
ber of sub-samples were used for comparison; 200 sub- 
samples are used in flgure[4]). The spectrum is peaked at 
a few frequencies, demonstrating how the wave pattern 
has a number of preferred wave- front velocities. This is 
not obvious from studying spatiotemporal diagrams such 
as Fig. [TJl. In the high frequency regime the power dies off 
in a manner consistent with a power law with exponent 
-2, consistent with the wave- fronts undergoing Brownian 
fluctuations. 



IV. CONCLUSIONS. 

We have studied granular pipe fiow by means of a one- 
dimensional two-parameter model. The two parameters 
are the coefficients of restitution for collisions between 
grains and collisions between grains and walls. The very 
simple collision rules are contrasted by the amount of 
structure exhibited by the model, and formation of den- 
sity waves varying greatly in magnitude and qualitative 
behavior is observed. We find a criterion for the forma- 
tion of density waves: that the dissipation from collisions 
with walls be greater than that from grain-grain colli- 
sion. Under this criterion our model predicts that den- 
sity waves can form also in the absence of any interstitial 
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gas. 

Contrary to intuition, the flow rate is largest when 
density waves are large, slow and far between. This indi- 
cates that in some circumstances, the flow rate in gravity- 
driven granular pipe flow can be increased by softening or 
roughening the pipe walls. For example, with soft grains 



described by /i = 0.3, a 'rough' pipe with v = 0.5 gives 
a flow rate two to three times faster than a 'smoother' 
pipe for which v = 0.7. 

We have benefited from discussions with J. P. Hulin 
and the Granular Media Group at FAST laboratory, Or- 
say. 
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